A new stability results for the backward heat 

equation 

Alain Pham Ngoc Dinh 

MAPMO, UMR CNRS 6628, Orleans University, France. 

Dang Due Trong 

Department of Mathematics ,HoChiMinh City National University, 
227 Nguyen Van Cu, Q. 5, HoChiMinh City, VietNam. 

Pham Hoang Quan 

Department of Mathematics, Sai Gon University, Ho Chi Minh city, Viet Nam 

Nguyen Huy Tuan 

Department of Mathematics and Informatics, Ton Due Thang University, 
98, Ngo Tat To, Binh Thanh district, Ho Chi Minh city , Viet Nam 

November 13, 2009 
Abstract 

In this paper, we regularize the nonlinear inverse time heat problem in the unbounded 
region by Fourier method. Some new convergence rates are obtained. Meanwhile, some 
quite sharp error estimates between the approximate solution and exact solution are pro- 
vided. Especially, the optimal convergence of the approximate solution at t = is also 
proved. This work extends to many earlier results in [3 [HI [101 [151 El OS 1211 124"] . 
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1 Introduction 



Transient heat conduction phenomena are generally governed by the parabolic heat con- 
duction equation and if the initial temperature distribution and boundary conditions are 
specified, then this, in general, leads to a well posed problem which may easily be solved 
numerically by using various methods. 

However, in many practical situations when dealing with a heat conducting body it 
is not always possible to specify the boundary conditions or the initial temperature. For 
example, in practice, one may have to investigate the temperature distribution and the heat 
flux history from the known data at a particular time. In other words, it may be possible to 
specify the temperature distribution at a particular time, say t — T > 0, and from this data 
the question arises as to whether the temperature distribution at any earlier time t < T 
can be retrieved. This is usually referred to as the backward heat conduction problem 
(BHCP), or the final boundary value problem. The systematic study of the backward heat 
conduction problem is of rather recent origin, although isolated considerations have been 
given to such problems for several hundred years. In general no solution which satisfies 
the heat conduction equation, the final data and the boundary conditions exists. Further, 
even if a solution did exist, it would not be continuously dependent on the boundary and 
the final data, see Payne [2]. Thus the BHCP is an example of an ill-posed problem that 
is impossible to solve using classical numerical methods and requires special techniques to 
be employed, see Hadamard [3] . Conditions for which the BHCP becomes well-posed have 
been investigated by Miranker [4] and Cannon and Douglas [5]. These studies introduced 
additional hypotheses which restrict the class of functions to which the solution must 
belong, and which are seldom satisfied. Therefore numerical methods of solution appear 
more useful. Thus regularization techniques, see for example Cannon [6] and Han et al. [7], 
have been developed for solving the BHCP. Different methods, based on a perturbation of 
the original parabolic heat equation were proposed by Lattes and Lions [8] and Lesnic et 
al. [9]. These methods consist of replacing the operator by a perturbed higher order one 
that holds better invertibility features. 

In this paper, we consider the following problem 

u t -u xx = 0, (x,t) G R x (0,T), 
u(x,T) = (p(x), x e R, 

where T > be a given and <p(x) are given. We want to retrieve the temperature dis- 
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tribution u(x,t) for < t < T. Of course, since the data </?(.) is based on (physical) 
observations, there will be measurement errors, and we would actually have as data some 
function (p t 6 L 2 (R) , for which \\(p — (p e \\ < e, where |.|| denotes the L 2 -norm, the constant 
e > represents a bound on the measurement error. That is to say, practically, we need to 
consider the following problem, 

u t -u xx = 0, (x,t) e Rx (0,T), 

u(x,T) = <fe(x), x G R, 

Notice the reader that the problem (1) is investigated in some recent papers of ChuLiFu 
0[8j[2l] and of other authors such as Lien[2T].Murniz [H], et al . To the authors knowledge, 
so far there are many papers on the backward heat equation, but theoretically the error 
estimates of most regularization methods in the literature are Holder type, i.e., 

\\u(.,t) -v e {.,t)\\ < Ce k , k>0. (3) 

where C is the constant depend on u, k is a constant is not depend on t, u . As we know, 
e k converges to zero more quickly than the logarithmic term. So, the major object of this 
paper is to provide new regularization method to established the Holder estimates such as 
([3]). We give a new approximation problem and investigate the error estimate between the 
regularization solution and the exact one. 

The remainder of the paper is divided into three sections. In Section 2, we establish the 
approximated problem and show that it is well posed. Then, we also estimate the error 
between an exact solution u of Problem (1) and the approximation solution v e (.,t) with 
the Holder type. Finally, a numerical experiment will be given in Section 3. 



2 Regularization and error estimates. 

Let 

i r+oo 

<p(0 = -= / ^x)e'^dx 

V Z7T J-r 



oo 

be the Fourier transform of the function (p e L 2 (R). We assume u(x,t) is the unique 
solution of ([1]). Using the Fourier transform technique to problem ([1]) with respect to the 
variable x, we can get the Fourier transform t) of the exact solution u(x, t) of problem 

'^(£,t) = (^)V£,*), 



The solution to equation (jlj) is given by 

u(Z,t) = e<?-*>?0(t). (5) 

or equivalently, 



1 



+00 



u(x,t) = —= / eV-^m^dt- ( 6 ) 

V ATI J 

—00 

Since t < T, we know from ([6]) that, when £ becomes large, exp (T — t)£ 2 increases rather 
quickly. Thus for u(£,t) G L 2 (R) with respect to £, the exact data function must 
decay rapidly as |£| — > 00. Small errors in high-frequency components can blow up and 
completely destroy the solution for < t < T.As for the measured data <p e {x), its Fourier 
transform <p e (0 is merely in L 2 (R). 

Let p > 1 be a constant number. We approximated problem (3) by pertubing the 
Fourier transform of final value <p as follows 



Of 



-T£ 2 



(7) 



The formal solution of (J7|) is also easily seen to be 



-t(. 



> 



*) = peiP-DTe+e-Te ^)' (8) 



or 
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we(x,t) = -±= J p e{p jj + e _ T p me* x d£, (9) 



where ft is a positive number such that ft > 0. 
Let f e (.,£) be the approximated solution given by 



1 f e -t « 

Ve(X,t) = 7^ J fte(p-^ + e -Te ^ e (10) 



Note that if ft is chosen small, then for small |£|, ^ e(p _ 1) e T g2 +e _ Te 2 i n @ is close to e^ T 
in (e. 
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For £, x, (3 > 0, < a < b, we prove the following inequality 

</rf. (n) 



e< 2 



1 + /3e 6 « 2 
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Thus, we have 



"■piiT^ " RF^ (12) 



e < 2 e < 2 



l + /3e fc € 2 " (l + /5e^ 2 )f(l + /3e b « 2 ) 1 -f 
e a 

< 



< /rf. 

For < t < s < T, denote 

exp{-^ 2 } 



B(Z,s,t) 



pexp{(p - 1)T£ 2 } + exp{-T£ 2 } 
exp {( s _t_T)e 2 } 



/?exp{(p- 1)T£ 2 } + exp{-T£ 2 }' 
Using the inequality ffTTT) . we obtain 



= TT^ 2 ^ (13) 
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B(t,s,t) = YV^-^' (14) 

Theorem 1. Let tp G L 2 (R). Then unique solution of the problem depends on the 
final value (p,i.e, ifw,v are the solution of the problem §6§ corresponding to the final values 
tp and u, then 

\\ w (.,t)-v(.,t)\\ Kp&y-ti 

Proof of Theorem 1. 

Let w and v be two solution of the problem (jSj) corresponding to the final values (p and (p. 
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Using Parseval inequality, we have 



\w(.,t)-v(., 



< 



\\w(.,t) - v(.,t)f 



fp(0 - 0(0 



Hence 



2t-2T ,, * ,, n 

< p— ||£-0|| 2 

2t-2T n 



Remark 1. 

In |H [19], the stability of magnitude is e~ . In [7J [151 120], the stability estimate is of order 



— -i 

e T 



>— -l 



T 



In our paper, we give a better estimation of the stability order, which is C(3t ( 
It is easy to see that the order of the error, introduced by small changes in the final value 
g, is less than the order given in |15l [20]. This is among of the advantages of our method. 
Theorem 2. Let w e and v e defined by (Q and ffTUl) . Then one has 



\w e {.,t)-v e (.,t)\\ <(3^e. 



Proof of Theorem 2. Using Theorem 1, we get 



\w e (.,t) - V t (.,t)\\ <(3vT \\ip e -ip\\ <(3pT 6. 



(15) 



(16) 



Theorem 3. Let u be the exact solution of problem such that ||ii(.,0)|| < E\. Then one 
has 



\\u(.,t)-v e (.,t)\\ < er(E 1 + l). 
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Proof of Theorem 2. Using (fill) , we have 

+00 

\\u(.,t)-w e (.,t)\\ 2 = [ \u(^t)-wp(Z,t)\ 2 d£ 



-00 
+00 



-00 
+00 



-00 
+00 



(1 + (3eP T Z 2 



,{pT-t)e 



Therefore 



1 + f3eP T t 

t — IT I 

< 2 [3^ J |£(£,0)f^ 
—00 

\u(.,t)-w e (.,t)\\ <PWE X . 



From (3 = e p and using the inequality, we obtain 



\u(.,t) - v e (. 



< \\u(.,t)-w e (.,t)\\ 

< 6?(E X + l). 



\w e (.,t) - v 6 (.,t)\\ 



Remark 2. Notice that the convergence estimate in Theorem 1 does not give any useful 
information on the continuous dependence of the solution at t — 0. This is common in the 
theory of ill-posed problems, if we do not have additional conditions on the smoothness of 
the solution. To retain the continuous dependence of the solution at t — 0, instead of (2.5), 
one has to introduce a stronger a priori assumption. 

We denote ||.||fc be the norm in Sobolev space H k (R), k > defined by 



IK,o)||, 



+00 



i + £TK£,o)|X 
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Theorem 3. Let u be the exact solution of problem such that ||w(.,0)||2 < E 2 . Let (5 = e , 
then one has 

VT t t-T+pT 

\\ u ^t)-w e (.,t)\\<f 7T -eWE 2 + e^-. (17) 
Proof of Theorem 3. 



\u(.,t)-w e {.,t)\\' 



p e (pT-t)e 



— oo 
+oo 

— oo 



2 +oo 



< i £lf& ] IK,o)||? 



ln( 



vT t 
ln (^) 



From p = e and using the inequality, we obtain 



\u{.,t)-v e {.,t)\\ < \\u{.,t)-w £ {.,t)\\ + \\w e {.,t)-v e {.,t)\\ 

pT _t_ t-T+pT 

ln(i) 



Remark 3. 



1. In t = 0, the error (1T7D becomes 



K.,0)- Ue (.,0 < J^-^ + eT. (18) 

It follows from p > 1 that the right hand side of (fTSI) converges to zero when e — > 0. This 
error is the same order as Theorem 3.1 in paper [8] (see page 566) and Theorem 2.1 in 
paper [9]. 

2. Since ffTTl) . the first term of the right hand side of (fTTj) is the logarithmic form, and 
the second term is a power, so the order of ffTTj) is also logarithmic order. It is the same 
order as some results which is mentioned in Remark 2. This often occurs in the boundary 
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error estimate for ill-posed problems. To retain the Holder order in [0, T], we introduce a 
different priori assumption. 
Theorem 4. 

Assume that there exist a positive number 7 G (0,pT) such that 

+00 



J e 2 ^u(Z,0)dt<El 



Let (3 = e and h = min{7, (p — 1)T}. Then, one has 



t + h 



\u(.,t)-w e (.,t)\\ <e^(E 3 + l). 



Proof of Theorem 4. 

From , we get 



\\u(.,t) -w e (., 



-00 
+00 



-00 
+00 



u^,t)-w^,t)\ 2 d^ 

!3 e vTe e (T-t)e 



-00 
+00 



[1 + f3eP T S 
Pe^ 2 



■m 



1 + (3eP T « 
1 + (3eP T ? 



d£ 



t+j-pT 



Using the inequality ( ITTi) . we have 

1 + f3eP T Z 2 
Combining ([20]) and (I2TI) . we get 

+00 

||u(.,f) - We (.,t)|| 2 < ffr J e 2 ^ 2 u(Z,0)dt < f^El 

—00 

From (3 = e and using the inequality, we obtain 

\\u(.,t)-v £ (.,t)\\ < \\u(.,t)-w e (.,t)\\ + \\w e (.,t)-v e (.,t)\\ 

t+7 t-T 
t + h 



(19) 



(20) 



(21) 
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Remark 4. 

1. In t = 0, the error (1191) becomes 

\\u(.,t)-w e (.,t)\\<eW(E 3 + l). (22) 

2. Suppose that E e = \\v € — u\\ be the error of the exact solution and the approximate 
solution. In most of results concerning the backward heat, then optimal error between is 
of the logarithmic form. It means that 

E t < C ^ln - 

where m > 0. The error order of logarithmic form is investigated in many recent papers, 
such as HE1 El [191 EH [22]. 

To illustrate this, we can enumerate some more recent papers considering the errors of 
logarithmic order. 

Let u and v e be exact solution and the approximated solution respectively. 
In [20[ [22] , the error is given the form 

\\v%,t)-u(.,t)\\<-^^. 

In recently, Feng Xiao-Li and coauthors [6] gave the error estimates as follows 

p 

" «(.)« < ^ + *iax{ [^rj 2 , a*}. 

From the discussed error, we can see that most of recent regularization methods established 
the logarithimic stability. The convergence rates in here is very slowly. 
In Tautenhahn and Schroter [T5] , the authors proved that the best possible worst case error 
for identifying u(t) is given by E t = De^ . 

In 2007, Zhi Qian et al. [8] (See Remark 3.6, p.570)gave the error estimation in t = for 
(2) as follows 

H.,0-v(.,0|| < - * . , s +max{l,T}/3*E. 
(ln(l/e))5 

where (3 = — -, — - — 

In (^(ln ir 

In [9], ChuLiFu and his coauthors esrablished the logarithimic order of the form 



\u(.,t) - u s>l 



, t i, E\ 2T ( ln E 

< E 1 -t In— 1 + 



& 



f lnf + ln(lnf )sr 
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Comparing (|22|) with the discussed error, we can see (122]) is the optimal error. 
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